Molecular Cell

How embryonic cells generate large clones of cells in the adult represents a fundamental question in biology. Here, using melanocyte stem cells (McSCs) in the zebrafish as a model, we explore the function of the master melanocyte transcription factor (MITF) in safeguarding McSCs in embryonic development and their potential to pigment large clones in the adult. MITF is well known is for its role in the specification of melanoblasts from the neural crest (NC) and their differentiation into melanocytes, yet little is known about how this activity shapes the stem cell lineages. Here, we use live imaging coupled with single-cell transcriptomics and lineage tracing to show that MITF (mitfa in zebrafish) protects the melanocyte stem cell (McSC) fate in zebrafish. Utilizing a temperature sensitive mitfavc7 mutant, we show loss of Mitfa leads to a surprising premature and aberrant expansion of McSC progeny at the niche during embryogenesis, coupled with novel emergent transcriptional cell states. Linage tracing of McSCs from the embryonic to juvenile stages reveals Mitfa activity is subsequently required in regeneration by Schwann cell-like and melanocyte stem cell progenitors that serve as a reservoir for fast-responding pigment progenitors. Thus, the impact of Mitfa loss on the melanocyte lineage is cell-state and stage-specific. The emergent cell states upon mitfa loss may have important implications for our understanding the loss of MITF activity in human genetic disease and melanoma.

Drug-tolerant persisters (DTPs) represent a major obstacle to durable responses in targeted cancer therapy. DTPs are commonly described as distinct single-cell states that survive drug treatment via reversible, non-genetic mechanisms and drive tumor recurrence. Recent work demonstrates that multiple DTPs can coexist, reflecting diversity in lineage, signaling programs, or stress responses. However, each DTP is still generally viewed as a uniform cellular phenotype. Building on our prior work describing a population-level DTP termed "idling" [Paudel et al., Biophys. J. (2018) 114, 1499-1511], here we present evidence supporting a fundamentally different view: that DTPs are not single-cell states, but rather heterogeneous populations composed of multiple sub-states with distinct division and death rates that balance to produce near-zero net population growth. Using single-cell transcriptomics and lineage barcoding, we identify multiple phenotypic states within idling DTP populations, with reduced heterogeneity compared to untreated populations, and find that idling DTP cells emerge from nearly all lineages. Transcriptomic and functional analyses further reveal altered ion-channel activity in idling DTPs, which we confirm experimentally. Moreover, drug-response assays reveal increased susceptibility of idling DTPs to ferroptosis, a non-apoptotic form of regulated cell death, indicating the emergence of vulnerabilities associated with drug tolerance. Altogether, our results support a population-level view of tumor drug tolerance in which DTPs comprise stable collections of phenotypic states, shaped by treatment-defined phenotypic landscapes, which are potentially vulnerable to subsequent interventions. This perspective implies that eradicating DTPs will require a fundamental shift away from cell-type-centric strategies toward sequential treatments that progressively reduce phenotypic heterogeneity by modulating the molecular and cellular processes that establish the DTP landscape, an approach previously termed "targeted landscaping."

Ewing sarcoma (EwS) is an aggressive, human-exclusive tumor typically driven by the EWS::FLI1 fusion protein. To assess whether the neomorphic functions of EWS::FLI1 are fundamentally dependent on evolutionarily recent cofactors such as ETS transcription factors (ETS-TFs), Plycomb group (PcG) proteins, CBP/p300, or specific subunits of the BAF complex, we expressed EWS::FLI1 in the model organism Saccharomyces cerevisiae. This minimal system was chosen because several key EWS::FLI 's cofactors possess greatly reduced sequence homology (e.g., BAF) or are lacking altogether (e.g., ETS-TFs, PcG, or CBP/p300). We used co-IP/MS to map the yeast interactome, Chip-Seq to identify gDNA binding sequences, RNA-Seq for global gene expression, and engineered reporters to test conversion of (GGAA) tandem repeats (GGAASat) into neoenhancers. We found that the yeast EWS::FLI1 interactome was more limited and qualitatively distinct from its human counterpart, sharing core machinery (e.g. RNA Polymerase II, FACT) but lacking the BAF/SWI-SNF and spliceosome complexes, and showing strong enrichment for the SAGA chromatin remodeling complex. We also found that EWS::FLI1 binds to hundreds of sites in the yeast genome with a clear preference for putative ETS-TF consensus sequences and (CA) dinucleotide repeats. Yet, EWS::FLI1 expressing cells presented only minimal transcriptional dysregulation, a stark contrast to the extensive changes observed in humans and Drosophila cells. Finally, we found that EWS::FLI1 successfully converted silent GGAASat sequences into active enhancers in yeast. This remarkable result occurs despite the absence of homologs for key human activators, such as CBP/p300, strongly suggesting that EWS::FLI1 can mobilize functionally related, non-homologous pathways to establish neoenhancers at GGAASat sites. Altogether, our results indicate that EWS::FLI1's core ability to drive GGAASat-dependent gene expression is a conserved, ancient property, while GGAASat-independent extensive transcriptome reprogramming is dependent on co-factors and pathways specific to animal cells.

Neuronal migration is a vital process that positions billions of neurons to create a functional brain. To navigate the constrained microenvironments within the cortex, precise control over the nuclear mechanics in migrating neurons is indispensable. Here, we show that Lamin B1 (LB1) regulates neuronal migration by modulating nuclear deformability. Excess LB1 in neurons halted migration without altering laminar identity or overall gene expressions in vivo, while in vitro, it elevated nuclear stiffness and impaired neuronal motility in confined spaces. Moreover, mispositioned neurons resulted in electrophysiological defects in the brain. Computational modeling predicted a temporal relationship between nuclear deformation and enhanced migration velocity, which was validated experimentally through live imaging. Notably, cerebral organoid assays using iPS cells established from patients with LMNB1 duplication exhibited impaired neuronal migration in a human model. Collectively, these findings demonstrate that LB1 is a critical regulator of nuclear mechanics, ensuring the accurate spatiotemporal positioning of neurons.

Patterns of locus coeruleus (LC) activity and norepinephrine (NE) release during non-rapid-eye-movement (NREM) sleep suggest a critical role for the LC-NE system in offline modulation of forebrain circuits. NE transmission promotes synaptic plasticity and is required for memory consolidation, but the field has only begun to uncover how LC activity contributes to coordinated forebrain network dynamics. Hippocampal ripples, a hallmark of memory replay, are temporally coupled with thalamocortical oscillations; however, the circuit mechanisms underlying systems-level consolidation across larger brain networks remain incompletely understood. Here, using multi-site electrophysiology, we examined LC firing in relation to hippocampal ripples in freely behaving rats. LC activity and ripple occurrence were state-dependent and inversely related: heightened arousal was associated with increased LC firing and reduced ripple rates. At finer timescales, LC spiking decreased {approx}1-2 seconds before ripple onset, with the strongest modulation during awake ripples but minimal change during ripple- spindle coupling. These findings reveal state-dependent dynamics of LC-hippocampal interactions, positioning the LC as a key component of a cortical-subcortical network supporting systems-level memory consolidation.

Circulating tumor cells (CTCs), and especially CTC-clusters, are linked to poor prognosis and may reveal mechanisms of metastasis and treatment resistance. Therefore, developing unbiased methods for the functional characterization of CTCs in liquid biopsies is an urgent need. Here, we present an evaluation of multiplex imaging mass cytometry (IMC) to analyze CTCs in mice with human xenograft tumors. In a single-step process, IMC uses metal-labeled antibodies to simultaneously detect a large number of proteins/modifications within minimally manipulated small volumes of blood from the tail vein or heart. We used breast cancer cell lines and a patient-derived xenograft (PDX) to assess antibodies for cross-species interpretation. Along with manual verification, HALO-AI-based cell segmentation was used to identify CTCs and quantify markers. Despite some limitations regarding human-specificity, this technology can be used to investigate the effect of genetic and pharmacological interventions on the properties of single and cluster CTCs in tumor-bearing mice.

De novo mutations (DNMs) arising in the parental germline are a major cause of severe developmental disorders. While most DNMs originate in the paternal germline, it remains unclear whether fathers of affected children carry a systematically altered burden of transmissible germline risk, or whether disease largely reflects stochastic outcomes of shared population-wide mutational processes. Here, we combined whole-genome sequencing of 168 parent-child trios with ultra-accurate duplex sequencing of paternal sperm to directly relate transmitted DNMs to the broader mutational and selective landscape of the male germline. In 127 fathers, sperm mutation burden and mutational spectra were indistinguishable from population reference cohorts. Positive selection metrics were likewise concordant, with a global dN/dS of 1.56 (95% CI 1.45-1.67) compared to 1.44 (95% CI 1.17-1.77) in controls and 28 of 32 significantly selected genes overlapping with prior findings. Six fathers harboured a pathogenic early mosaic variant detectable in sperm at allele fractions that ranged from 0.7% to 14.8%. Although these variants generated substantial individual-level risk outliers, they accounted for only [~]11% of the aggregated exome pathogenic burden across the cohort. The remaining burden was distributed across low-VAF mutations, including positively selected driver variants and other rare mutations accumulating with paternal age. Together, these results show that transmissible de novo disease risk is governed primarily by universal germline mutational and selective processes, while early developmental mosaicism produces uncommon but clinically meaningful deviations. This integrated view clarifies how mutation timing, age-associated accumulation and germline selection jointly shape inheritance risk.

Pathogenic variants in TP53, the key tumour-suppressor gene underlying Li-Fraumeni syndrome (LFS), are among the best-established causes of inherited cancer predisposition. However, large-scale sequencing has revealed that many apparently pathogenic TP53 variants detected in blood are the result of somatic clonal expansions, complicating risk interpretation. Using blood-derived whole-exome data from 469,391 UK Biobank participants, we combined variant allele fraction (VAF) with haplotype-sharing analysis to distinguish germline and somatic TP53 variants. Germline variants were concentrated at sites linked to partial loss of p53 function and lower disease penetrance, whereas classic LFS alleles appeared almost entirely somatic. High-VAF carriers of classic LFS alleles conferred markedly increased risk of haematological malignancy but not solid tumours, consistent with large TP53-mutant clonal expansions. The prevalence of somatic clonal expansion also correlated with missense variant pathogenicity, suggesting that somatic activity provides an informative in vivo proxy for functional impact. These results provide new insights into TP53-associated cancer risk at the population level, demonstrate that somatic rather than germline risk predominates in middle-aged healthy adults and provide a scalable framework for variant classification in large-scale population genomics.

Systemic inflammation is implicated in various diseases, yet its upstream determinants remain poorly examined. We conducted a large scale two-sample Mendelian randomisation (MR) study to systematically evaluate the potential causal effects of 3,213 molecular (metabolomic, proteomic), physiological and disease traits on circulating interleukin-6 (IL-6) and C-reactive protein (CRP) levels. Genetic instruments were derived from genome wide association studies and analysed using inverse variance weighted (IVW), weighted median, and MR-Egger methods with multiple testing correction. Bidirectional MR was performed to assess reverse causation. After Bonferroni correction, evidence of potential causal effects was observed for 72 traits on CRP and 9 traits on IL-6. CRP was predominantly influenced by metabolomic traits, especially lipid and fatty acid measures. Genetically proxied adiposity (body mass index and obesity), triglyceride rich lipoproteins, glycoprotein acetyls (GlycA), and apolipoprotein E increased CRP levels, whereas HDL-related cholesterols, polyunsaturated fatty acids, and glutamine decreased CRP. Most associations were consistent across MR methods, supporting the robustness of these results. As expected, IL-6 had a large effect on CRP. IL-6 was influenced by primarily adiposity and HDL-related lipid measures, with generally smaller effect sizes and limited support across sensitivity analyses. Bidirectional analyses indicated little evidence that CRP directly drives metabolic traits when restricting to cis-acting instruments, whereas genetically proxied IL-6 signalling showed consistent downstream effects on HDL particle concentration and composition. Adiposity is a shared upstream determinant of both inflammatory biomarkers, with stronger and broader effects on CRP. These findings suggest that CRP acts as an integrated downstream readout of systemic inflammatory burden, whereas IL-6 reflects a more tightly regulated and context-dependent process. Our work clarifies traits that may causally influence systemic inflammation and highlights biological pathways linking inflammation to cardiometabolic and inflammatory diseases. By mapping upstream determinants of IL-6 and CRP, we also provide a resource to prioritise key drivers for mechanistic study and therapeutic targeting.

Quantitative measures for tracking functional health have generally been lacking. Intrinsic capacity (IC) has been proposed as an appropriate measure, but its metrics have been derived in small datasets and sparse longitudinal data. Using harmonized measures of cognition, locomotion, sensory function, vitality, and psychological well-being from 501,615 UK Biobank participants and followed for a median of 15.5 years, we derived domain-specific and composite IC scores. We examined associations with incident disease, cause-specific mortality, multimorbidity, lifestyle and socioeconomic factors, and multi-omic profiles from Olink proteomics, NMR metabolomics, clinical biochemistry, and blood-cell traits. We found that composite IC declined non-linearly with age, and within-person decline was steeper than the cross-sectional age measures. Participants with greater baseline morbidity, those who subsequently developed incident disease, and those who died earlier in follow-up showed lower IC trajectories across adulthood. The IC domains were only modestly correlated with one another, supporting multidimensionality, yet higher overall IC was associated with lower risk of most diseases examined. The dominant IC domain varied by endpoint, with cognition informative for dementia, sensory function for hearing loss, psychological capacity for depression, locomotion for osteoarthritis, and vitality for cardiometabolic outcomes. IC was also associated cross-sectionally with physical activity, insomnia, smoking, medication burden, and socioeconomic disadvantage. More proteins were found predictive for vitality, and enrichment converged on immune/inflammatory and metabolic pathways. Blood-based surrogates recapitulated part of the phenotypic signal, particularly for vitality. Overall, this IC framework captures longitudinal health trajectories and broad disease vulnerability in a large middle- to older-aged cohort and supports IC as a clinically meaningful, multidomain phenotype of aging and identifies blood-based correlates that may facilitate at-scale future monitoring of aging-related function declines.

Dementia affects tens of millions of people worldwide, yet disease-modifying treatments remain strikingly limited. Although the recombinant zoster vaccine Shingrix has been associated with reduced dementia incidence, its potential influence on individuals already living with dementia is unknown. Here, we followed a propensity-score matched cohort of 68,960 US dementia patients using a nationwide electronic health record network, comparing Shingrix recipients within two years of diagnosis to recipients of any other vaccine. Shingrix was associated with substantially reduced all-cause mortality across the first three years of follow-up (hazard ratios 0.74, 0.88, and 0.89; P[≤]0.006), robust across multiple sensitivity analyses. Furthermore, within-individual subgroup analyses of repeated Mini-Mental State Examinations conducted 3-6 years apart revealed significantly divergent cognitive decline rates across groups (time-by-group interaction P=0.002). Interval vaccination was associated with more stable cognition, contrasting with steeper declines in unvaccinated individuals. These findings support prospective evaluation of recombinant zoster vaccination as a potential strategy to improve outcomes in patients with established dementia.

The transcription coregulator OCA-B promotes CD4+ T cell memory recall responses and autoimmunity. OCA-B T cell deletion prevents spontaneous type-1 diabetes (T1D) onset in non-obese diabetic (NOD) mice and blunts T1D in a subset of more aggressive models. However, the role of OCA-B in diabetes induced by treatment with immune checkpoint inhibitors (ICIs), and the role of OCA-B in the control of tumors with and without ICI treatment, has not been studied. Here we show that islet and pancreatic lymph node T cells from T1D individuals express measurable POU2AF1 mRNA. Deletion of OCA-B in T cells fully insulates 8-week-old non-obese diabetic (NOD) mice against ICI-induced diabetes and partially protects 12-week-old mice. Salivary and lacrimal gland infiltration and inflammation were also reduced. Protection was associated with a block in the differentiation of progenitor exhausted CD8+ T cells (TPEX) into terminally exhausted CD8+ T cells (TEX). We show that OCA-B T cell loss preserves anti-tumor immune responses following PD-1 blockade in different tumors and mouse strains. These findings point to a potential therapeutic window in which pharmaceuticals targeting OCA-B could be used to block the emergence of both spontaneous and ICI-induced autoimmunity while sparing anti-tumor immunity. We develop first-in-class small molecule inhibitors of Oct1/OCA-B transcription complexes and show that administration into NOD mice also blocks diabetes emergence following PD-1 blockade. These results identify OCA-B as a promising therapeutic target for the prevention of autoimmunity and immune-related adverse events (irAEs).

Public RNA-seq sample sets can refine per tumor diagnosis and risk, but heterogeneous biology and analytic drift often obscure structure. Dynamic Quantum Clustering (DQC), an unsupervised geometry-preserving method requiring no clinical labels or preset cluster counts, addresses both challenges. Applied to RNAseq from 692 TCGA gliomas (524 low-grade gliomas (LGG), 168 glioblastomas (GBM); 20,057 protein coding genes), DQC produced two dominant clusters with 90.9% post hoc diagnostic concordance and clear survival time separation. Filtering genes by inter-cluster mean differences yielded a 554 gene subset that improved accuracy to 97.3%. Rank ordering these genes identified ~90 genes that, under DQC, produced three LGG-pure subclusters with ordered, but different survival outcomes and one GBM-rich cluster (PPV 97.1%)--the RNA-based clustering without clinical information thereby inherently reveals molecular groupings which mirror critically important clinical features. Comparing these clusters defined four nonoverlapping gene modules and assigned four BioCoords per tumor. DQC with Biocoords recapitulated the LGG-to-GBM continuum with a mesenchymal/invasion-extracellular matrix axis exhibiting a monotonic survival gradient, illustrating how geometry-aware unsupervised learning can translate bench and computational discovery into meaningful biology-based patient stratification and prognosis.

The APOE {varepsilon}4 allele is the strongest genetic risk factor for Alzheimers Disease. However, its distribution across Indian populations is poorly characterized. We analyze APOE allele frequencies in 9,524 individuals from 83 distinct populations in the GenomeIndia dataset. {varepsilon}4 frequencies show large variation across populations within India, ranging from 2.7% to 36.1%, with a median of 11%. Tribal populations have higher {varepsilon}4 frequencies compared to non-tribal groups, while Tibeto-Burman populations have significantly lower frequencies. One tribal population from the northern coastal highlands has {varepsilon}4 frequency of 0.36, with 59% of individuals being carriers. {varepsilon}4 carrier status correlates significantly with lipid phenotypes including LDL, HDL, total cholesterol, and triglycerides. Collectively, these findings reveal exceptional genetic diversity in Alzheimers Disease risk across India and have important implications for population-specific screening strategies, genetic counseling, and precision medicine approaches to dementia prevention.

Children with acute lymphoblastic leukemia (ALL) are treated with multiagent chemotherapy that causes profound changes to the immune system. There are limited data on how disease and therapy impact antigen-specific immune memory, leading to inconsistent guidelines on best practices for revaccination of this population. Here, to inform vaccine guidance, we investigated whether immunity derived from routine childhood measles and varicella zoster virus (VZV) vaccines is maintained during and after therapy for childhood ALL. We report that antibodies against measles and VZV were significantly reduced in children with ALL (n=45) compared to healthy controls (n=13), particularly in older children in whom a longer time had passed since their most recent vaccine dose. However, the avidity of the measles and VZV-specific antibodies was indistinguishable between groups. Despite changes to the composition of the T cell compartment, both overall and antigen-specific T cell function were preserved in children with ALL. These data provide compelling evidence for revaccination of children following ALL treatment. Intact T cell responses suggest that post-treatment revaccination would be effective.

The shared genetic liability between cortical morphology and psychiatric disorders remains unclear. We aimed to identify whether the genetic loci shared between cortical morphology and six psychiatric disorders show regional or global effects. We identified substantial pairwise genetic overlaps of cortical surface area or thickness with psychiatric disorders; however, these loci lacked a uniform direction (~50% concordance). Cross-trait analyses revealed distinct architectures: internalizing disorders and schizophrenia/bipolar disorder shared more genetic loci with localized effects, whereas neurodevelopmental disorders shared fewer loci but more with widespread effects. We identified 17 genomic loci shared across all disorders, most of which had opposing directional effects across cortical regions. Only one locus (rs2431112) had region-specific and unidirectional effects (reduced primary visual and posterior cingulate surface area). This directional heterogeneity within and across pleiotropic loci reveals complex shared genetic architectures and likely limits the genetic predictive performance of brain morphology for psychiatric disorders.

Mitral valve prolapse (MVP) is the most common cause of primary mitral regurgitation and is associated with the development of malignant arrhythmias, often in the context of myocardial fibrosis. The genetic architecture of MVP, and whether there are genetic factors explaining why only some individuals with MVP have adverse outcomes, remains poorly understood. We performed a meta-analysis of genome-wide association studies (GWAS) for MVP encompassing 21,517 cases among a total sample size of over 2.2 million individuals. We discovered 89 genomic risk loci for MVP, of which 72 were novel findings. Prioritization of causal genes and pathways using epigenetic and transcriptomic data from mitral valve and extra-valvular tissues replicated known gene associations to MVP including those involved in TGF-{beta} signaling and extracellular matrix biology, but additionally emphasized a role in MVP for biological pathways relevant to cardiomyocyte biology. Accordingly, we identified several MVP risk loci with pleiotropy to cardiomyopathies, especially hypertrophic cardiomyopathy, and demonstrated a significant genetic correlation between MVP and hypertrophic cardiomyopathy. Finally, we interrogated snRNA-seq data in human papillary muscle tissue from two individuals with severe MVP, characterizing genes associated with both risk of papillary muscle fibrosis and MVP.

ELF4 is an ETS family transcription factor involved in immune regulation, and germline loss-of-function mutations in ELF4 have been known as deficiency in ELF4, X-linked (DEX). To date, ELF4-related disease has been exclusively associated with germline mutations. Here, we report a pediatric patient with recurrent mucocutaneous inflammation and periodic fever caused by a somatic truncating mutation in ELF4. By directly comparing ELF4-mutant and wild-type immune cells within the same individual using full-length single-cell RNA sequencing, we identified mutation-associated transcriptional alterations across multiple immune cell types. Pathway analyses revealed cell type-specific immune alterations, characterized by reduced antiviral and interferon-related signaling in NK cells and enhanced inflammatory pathways related to Th17 differentiation and inflammatory bowel disease in CD16 monocytes. This study expands the disease spectrum of ELF4 deficiency by identifying somatic truncation of ELF4 as a genetic mechanism underlying autoinflammatory diseases and biased immune programs.

Our understanding of human mucosal T cell clonotype distribution in health and disease has centered on immunodominant antigens. We performed single cell T cell receptor (TCR) and RNA sequencing as an untargeted approach to define distributions of T cell clonal groups in health and ulcerative colitis (UC) across 333,088 T cells in colon and peripheral blood. Healthy donor-specific TCR repertoires had limited blood-colon clonal sharing, which was highest in cytotoxic T effector memory (Tem) populations and lowest in regulatory T cells (Tregs), reflecting tissue-based compartmentalization. Within healthy colon, TCR repertoires showed high T cell clonal sharing independent of anatomic distance, associated with high intra-clonal phenotypic diversity. Colon cytotoxic and Th17 populations showed high dispersion across sites, while Tregs were compartmentalized. Clonal lineages dispersed across blood and colon upregulated trafficking markers, suggesting active movement between tissues, while those dispersed across colon sites upregulated residency markers, suggesting intra-colon repertoire sharing is mediated by long-term, slow moving clonal groups. In UC, Tregs were expanded across inflamed sites, and increased CD8 Tem clonal groups showed increased dispersion regardless of inflammation. These findings reveal principles of T cell clonal organization in the human colon during health and disease, identifying opposing patterns of clonal dispersion among Treg and Th17 clonal groups, high phenotypic diversity within dispersed clonal groups, and elevated cross-colon dispersion of CD8 Tem clonotypes in UC.

Targeted protein degradation (TPD) is a therapeutic strategy to remove disease-causing proteins by routing them to the ubiquitin-proteasome, autophagy, or lysosme machineries. For instance, proteolysis-targeting chimeras (PROTACs) are synthetic hetero-bifunctional small molecules that simultaneously bind the target and an E3 ubiquitin ligase to drive ubiquitination and degradation by the proteasome. Despite considerable success, designing such molecules is challenging and the number of currently addressable ubiquitin E3 ligases is limited. Here we demonstrate hetero-bifunctional de novo designed proteins as alternatives for TPD to access more targets and ligases. First, we develop a stable and highly adaptable helix-turn-helix scaffold for presenting different binding sites. Next, we use computational protein design to incorporate and embellish hot-spot- binding sites to target BCL-xL, plus short linear motifs (SLiMs) for KLHL20 ligase recruitment. The resulting mono- and bi-functionalised proteins bind the targets in vitro, and the latter degrade BCL-xL in cells leading to apoptosis.